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Abstract: We study the structure of a family of static, spherically symmetric space-times 
generated by an anisotropic fluid and governed by a particular type of f{R) theory. We 
find that for a range of parameters with physical interest, such solutions represent black 
holes with the central singularity replaced by a finite size wormhole. We show that time-l ik e 
geodesics and null geodesics with nonzero angular momentum never reach the wormhole 
throat due to an infinite potential barrier. Eor null radial geodesics, it takes an infinite 
affine time to reach the wormhole. This means that the resulting space-time is geodesically 
complete and, therefore, nonsingular despite the generic existence of curvature divergences 
at the wormhole throat. 
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1. Introduction 

Our understanding of the physical world is nowadays based on two basic pillars, namely, quantum 
theory and the notion of curved space-time. The enormous individual success of these approaches 
contrasts with the difficulties that arise when one attempts to combine them. This is particularly 
evident when one explores quantum properties of black holes. Despite the classical robustness of these 
objects, which rapidly decay after they are formed into stationary configurations characterized by their 
mass, charge, and angular momentum [1-3], Hawking discovered that they are affected by a quantum 
instability due to the very existence of an event horizon [4], with puzzling implications for the unitary 
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evolution of quantum states [5]. This low energy phenomenon is sometimes regarded as the tip of a larger 
ieeberg that dwells deeper down the blaek hole. The emitted Hawking quanta, whieh possess positive 
energy and eause a slow evaporation of the blaek hole, are eorrelated with negative energy modes that 
fall into the eentral singularity. The preeise fate of the infalling quanta is eompletely unknown due to the 
breakdown of the geometrie deseription as the singularity is approaehed. An improved understanding of 
the physieal proeesses taking plaee in the innermost regions of blaek holes eould thus help shed some 
new light on the evaporation proeess and the interaetion between geometry and the quantum. 

As a way to get a new perspeetive on the internal strueture of blaek holes, we initiated a researeh 
program based on the notion of eurved spaee-time but allowing the geometry to be slightly more general 
than the standard Riemannian framework of elassieal General Relativity. In this program we allow 
the affine eonneetion to be determined by a variational prineiple rather than being fixed a priory by 
the Christoff el symbols of the spaee-time metrie. This approaeh is supported, mainly, by the different 
mathematieal and physieal meanings that metrie and affine eonneetion have [6,7], and by the empirieal 
evidenee of the emergenee of this type of geometries in eondensed matter systems with mieroseopie 
defeets [8]. If the spaee-time eontinuum that we pereeive had a mieroseopie strueture, a lesson that 
eould be extraeted from eondensed matter systems is that metrie-affine geometries eould be the natural 
framework to formulate gravitational theories and their phenomenology [9] . 

By reeonsidering the roles typieally attributed to metrie and affine struetures in eertain extensions 
of General Relativity (GR), we have been able to show that blaek hole singularities ean be avoided in 
eertain elementary eonfigurations. To be more preeise, we have eonsidered eleetrovaeuum spaee-times 
in theories of gravity ineluding quadratie eurvature terms [10-15] as well as in others inspired by the 
Born-Infeld nonlinear theory of eleetrodynamies [16,17]. The novel and main result of our analysis is 
that spaee-time singularities ean be avoided without neeessarily removing eurvature divergenees [18]. 
The absenee of singularities in these spaee-times is due to the faet that the GR point-like eentral 
singularity is generieally replaeed by a finite-size wormhole strueture. We point out that sueh wormholes, 
as opposed to those typieally eonsidered in the framework of GR [19], do not violate any of the energy 
eonditions. This meehanism to remove spaee-time singularities is in sharp eontrast with the standard 
state-of-the-art of the field, where the widespread identifieation between eurvature divergenees and 
spaee-time singularities has developed multiple strategies foeused on bounding eurvature sealars as a 
way to obtain nonsingular eonfigurations [20-24], a program that earries its own limitations, sueh as the 
diffieulty to find theoretieally well-grounded theories able to support sueh nonsingular solutions. To the 
light of our results, this identifieation has to be reeonsidered, and further exploration of extensions of GR 
with independent metrie and affine struetures is needed. 

In this work we show that even in a simpler seenario, sueh as in the well-known f{R) theories of 
gravity (see [25] for a review), one ean find blaek hole eonfigurations whieh represent nonsingular 
spaee-times. Classieally, singular spaee-times are eharaeterized by the existenee of inextendible 
paths [26-28] (see also [29] and referenees therein). In partieular, if physieal observers (represented 
by time-like or null geodesies) start or terminate their paths at a finite value of their affine parameter, 
then the spaee-time is said to be singular. If those geodesies ean be extended to arbitrarily large values 
of their affine parameter in the past and in the future, then the spaee-time is nonsingular. 
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Like in the more sophistieated models eited above, in the f{R) case considered here the extendibility 
of geodesics is guaranteed by the existence of a robust wormhole geometry at the center of the 
black hole. As we will see, this wormhole is supported by an energy-momentum tensor with the 
form of an anisotropic fluid which satisfies all the energy conditions. Interestingly, we find that the 
resulting geometry near the wormhole is functionally identical (up to constant coefficients) to that 
found in configurations generated by an electric field governed by the Born-Infeld nonlinear theory 
of electrodynamics [15]. This fact allows us to conclude that the wormhole structure is a genuine and 
robust property of these theories, having very little sensitivity to the details of the particular matter source 
chosen. In the asymptotic region far away from the wormhole, however, the geometry is very dependent 
on the particular details of the matter source. As we shall show, this is due to the fact that it is possible 
to establish a direct correspondence between the anisotropic fluid used here and a certain family of 
nonlinear theories of electrodynamics, which reinforces the original motivation of [15] to explore f{R) 
theories coupled to that particular type of matter. 

2. Field Equations and Matter Source 

To fix ideas, the action we are interested in can be written as 

s = ^ j d'^Xyf^f{R) + ( 1 ) 

where is related to Newton’s constant (in GR, n? = SttG), g is the determinant of the space-time 
metric f{R) is a given function of the curvature scalar R = gf^'^Ri^i^{T), and is the matter action 
depending on g^i, and on the matter fields The key idea in our program is that we allow metric and 
affine connection T = T^^^ to be independent (Palatini or metric-affine approach), which means that both 
of them are governed by the variational principle. Using the fact that T appears in the action Equation (1) 
only through the Ricci tensor R^,y{T), independent variations of (1) with respect to g^^ and T^^ yield the 
field equations [25] 


f rR^u 2 



( 2 ) 

(3) 


where = —^^5Sm/Sgij.i, is the energy-momentum tensor of the matter. It is well known that the 
connection Equation (3) can be formally solved by introducing a conformally related metric = fug^i, 
such that this equation becomes V= 0, which implies that the connection is compatible 
with i.e., that it can be written as the Christoffel symbols of hence i?^y(r) = R^y{h). With 

this result, the remaining field Equation (2) for Palatini f{R) theories can be written in the following 
compact form: 


R 


V- (^) ~ p 

^ j 


R L 


^ 5 / + 


(4) 

Here we are denoting R^^{h) = R^a{h)h'^’^, with R^a{h) the Ricci tensor of the metric One must 
bear in mind that tracing over Equation (2) with the metric g^'^ one gets Rf^ — 2/ = r^T (with T the 
trace of the energy-momentum tensor), which is an algebraic relation leading io R = R{T). This implies 
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that all the objects on the right-hand-side of Equation (4) are functions of the matter. We also note that 
Equation (4) represents a system of second-order equations for (and also for because they are 
algebraically related). The second-order character of the field equations is a robust property of Palatini 
theories [25,30-34]. 

Given that the conformal factor that relates with is a function of the trace T of the matter 
sources, for traceless sources such as a standard Maxwell electromagnetic field the metric field equations 
boil down exactly to those of GR for arbitrary f{R) theory. To avoid this degeneracy and explore new 
dynamics in black hole configurations, in [15] we considered matter sources represented by nonlinear 
theories of electrodynamics, where the Maxwell action is replaced by a given function of the gauge field 
invariants. These theories yield, in general, a nonvanishing trace in four dimensions. Different new 
black hole configurations were found in that work, including one with unusual properties corresponding 
to the case of f{R) = R — \R^ coupled to the Born-Inf eld theory of electrodynamics [16]. At that time, 
we simply noticed that the central singularity had been shifted from the origin to a finite radius. With 
our subsequent analyses in other gravity theories, we can now revisit those results and interpret them 
under the light of a new family of black holes, namely, those with a finite-size wormhole replacing the 
point-like central singularity. We will see next that the same geometric properties found in the case of 
having a nonlinear electric field as matter source are also reproduced by an anisotropic fluid of the form 

= diag[-p,P^,P0,P^] 

= diag[-p, -p, ap, ap] (5) 

with a a constant, being a = 1 equivalent to the case of a Maxwell electric field. To satisfy the 
energy conditions, one requires 0 < a < 1. This fluid has recently been used in [35] to generate an 
interesting one-parameter family of wormholes in the framework of Bom-Inf eld gravity [36,37], and 
also considered in [38] to discuss generic conditions for generation of wormholes (let us also point out 
that such anisotropic fluids, which are believed to be relevant for the description of compact objects, can 
be modelled in terms of a perfect fluid, a standard electromagnetic field, and scalar field [39], which 
could have relevant implications for the description of spherically symmetric space-times). We will see 
that something similar occurs in the case of /(P) gravity, confirming in this way that the wormhole 
structure of the /(P) = P — AP^ model (and possibly others as well) is rather insensitive to the details 
of the matter sourcing it, thus supporting its robustness. 

3. Structure Equations 

To proceed, we consider a line element for of the form 

+ . ] . -f f‘^{x)dVP‘ (6) 

^ ^ A{x) ^ ’ 

where 'il){x) and A{x) are functions to be determined by solving the field Equation (4), which can be 
written as 

(2 - ^ ^ A 

/I \ 0 (I-f a/t^p) J j 


R,%h) 


(7) 
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where J and 0 are the identity and zero 2x2 matriees, respeetively. This strueture of the field equations 
implies that Rt — Rx^ = 0 , whieh allows to set 7p{x) —0 and P{x) —)■ without loss of generality, 
turning the auxiliary line element Equation ( 6 ) into a standard Sehwarzsehild-like form 

ds'^ = —A{x)dt^ + . dx^ + x^dVt^ ( 8 ) 

^ ^ A{x) 

The line element for the spaee-time metrie ean thus be written as 


ds^ = (—A{x)dt'^ + . dx‘^] + r'^{x)dQ? (9) 

Jr V Mx) ) 

where we have defined r'^{x) = x^ j f^. Using this line element, the eonservation equation of the fluid 
Equation (5) leads to p{x) = [35], where U is a eonstant with dimensions (in the a = 1 

ease, C would be proportional to the eleetrie eharge squared). With this result, we find that for the model 


f{R) = R-XR^ (10) 

we get R = 2(1 — a)n?‘p{x), and fn = 1 — (4A)k^( 1 — Q;)C'/r(a:)^+^", whieh ean be written 
as fR = 1 — by defining r = TcZ and = (4 A)k^(1 — a)C. Now, using the Ansatz 

A{x) = 1 — 2M{x)/x, from the equation Re^ we get 


8 ^ .X 

-T^x = 

' n 


(l-g 2^2+2°) 


( 11 ) 


It is now eonvenient to reseale x as x = rpx to obtain 


X = Zy/l - l/z2+2a (12) 

and dx = [(1 + aj— 1/z‘^+‘^°‘]dz, and to write M(x) as M(x) = Mq + Mq5iG{x), with 
= r^/ (SXMq). With all this. Equation (11) turns into 




(l + ^2+2c) 1 _ 1 

^ 2 a 1 — a 2 ^ 2 + 2 “ 

J R 


This equation ean be readily integrated, yielding 


(13) 


4 q ,_1 


G^{z) = 


z'^‘^+'^^/lA^^RG+r)(^2a‘^+a+2z^°‘+^-3) 8a(a2-l) 


2 ; 2 a + 2 _]^ 


4q;H-1 


2(a - l)(2a - 1) 


(14) 


where 2 F 1 {R,b,c;y) is a hypergeometrie funetion. The above solution is valid for a 7 ^ 1,1/2. As 
mentioned earlier, when a = 1, the fluid we are eonsidering degenerates into an eleetrie Maxwell field, 
thus yielding the Reissner-Nordstrom solution of GR (with fn = 1, = x^, and G{z) = —1/z). Eor 

a = 1 / 2 , the integration of G^ must be performed separately and the result is 


Gi 


— %^/zAAFIz^R log (z^R + y/z^ — 1 ) + 1 

7- - 


(15) 
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where 7 = | — is a eonstant neeessary to mateh the GR solution in the region 2 ; 3 > 1 (in GR 

Gi{z) = 21n2;). For eompleteness, we mention that in GR we have G^^{z) = (^i_a){i- 2 a) ^ whieh is 
reeovered from Equation (14) in the z ^ 1 limit. 

We note that the GR limit is reeovered very quiekly, just a few units away from z = 1, whieh implies 
that the number and loeation of horizons in these solutions is almost eoineident with those predieted by 
GR. In this sense, one ean find eonfigurations with two horizons, with one degenerate horizon, or without 
horizons, depending on the partieular eombination of parameters (a, ^ 2 ) chosen. 


4. Wormhole geometry. 


As found in other Palatini theories [10-12,17], the existence of a wormhole can be inferred from the 
minimum of the area function A = 47 rr^(a:) at a: = 0 (where z = 1). The explicit relation r = r(x) 
can be obtained for each a by inverting Equation (12). Eor several values of a the result is plotted in 
Eigure 1 . Near the wormhole throat, the metric component 


9tt 


1 + 6iG{z) 

hzfT 


/ Ir 


(16) 


with S 2 = Tc/ 2 Mo can be expanded as follows 


9tt 




47^a(»+l)^ir(|g±f) 
(8«2 2a-l)r(|^) 


852(1 






{2a + 2 ) 3/252 {z — 1)^^^ 


+ 


(3 — 6a)5i + 4(a — 1)52 1 


8(a2-l)52 (-2-1) 


+ 


(6a + 5) ((8a^ - 2a - 1) E (|g±|) - 4^a{a + l)5ir (|g|)) 
8V^^^ {V2{a + l)3/2(2a - l)(4a + 1)521 (g^)) 

(—244a^ + 4a + 171) 5i + 64 (2a^ + a — 3) 52 , , - 

-256(a-l)(a + lM,-+ 


(17) 



Figure 1. Representation of the radial function r{x) in units of the scale Tc for the cases a = 
1/10,1/2,4/5 (colors blue, green, and orange, respectively). The dashed lines represent the 
function |a:|, corresponding to the General Relativity (GR) case, which is quickly recovered 
for X > 2. As can be seen, the dependence on a is very weak. 
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This expansion is functionally identical to that found in [15] for this gravity model 
coupled to the nonlinear Born-Infeld theory of electrodynamics, with divergences of the form 
Qtt ~ ~ lY -\-h/{z — 1)^/^ + ..., which induces similar divergences in the Kretschmann scalar, 

K = with leading order K ~ — 1)^. 

The behavior far from the wormhole, however, does depend on the kind of matter source. In fact, in 
the Born-Infeld case of [15], the asymptotic geometry recovers the usual Reissner-Nordstrom solution, 
Qtt = — ill the notation of this paper, whereas in the case of the fluid considered here, 

the asymptotic behavior is 

^ ^ <52(l-«)(l-2«)z2a) 

which does depend on the specific fluid. This result can be understood as follows. Consider an 
electromagnetic field described by the action 

= ^ j d^x^/^ip{X) (19) 

where (p{X) is a given function of the electromagnetic field invariant X = where 

= d^Ay — dyA^ is the field strength tensor of the vector potential A^y. Different functions '^{X) 
specify different models of nonlinear electrodynamics. For electrostatic configurations, F^'' = E{r), the 
components of the energy-momentum tensor become 

= ^diag[<y9 -2Xipx,(p -2Xipx,(p, (f] (20) 

If we impose the equivalence of this energy-momentum tensor with that of the anisotropic fluid 
Equation (5), one finds the following function 


(^(X) = XP (21) 

where p = (1 -f a)/(2Q;), as the matter nonlinear electromagnetic Lagrangian. This corresponds to a 
class of models falling within the general analysis carried out in [40,41] for gravitating electromagnetic 
fields within GR. In addition, the properties of Equation (21) in different gravitational scenarios have 
been studied in detail in [42^6]. For our analysis, the relevant point here is that the asymptotic 
behavior of gravitational configurations sourced by the models Equation (21) is governed by the standard 
Schwarzschild mass term if 1 < p < 3/2, corresponding in our case tol/2<Q!<l [see Equation (18)], 
while those satisfying p > 3/2 (0 < a < 1/2) are asymptotically flat but not Schwarzschild-like, 
since they are governed by the charge term. In the case p = 1 (a = 1) we get a standard Maxwell 
field and thus the standard asymptotic Reissner-Nordstrom behaviour of GR is recovered. This is the 
case of Born-Infeld electrodynamics mentioned above. Nevertheless, these asymptotic behaviors have 
no influence on the features of the wormhole structure, since this arises due to the interplay between 
the f{R) gravity and the anisotropic fluid (nonlinear electromagnetic field) as the innermost region is 
approached for every a, which indicates that the wormhole geometry is a robust structure insensitive to 
the matter source that generates it. 

Eet us point out that the existence of a wormhole structure seems to be a robust prediction of Palatini 
gravities, since they are also present in quadratic [10-14] and Born-Infeld gravity theories [17] sourced 
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by (linear and non-linear) electromagnetic fields. The emergence of wormholes in scenarios with 
self-gravitating fields allows to interpret such solutions as geons in Wheeler’s sense [47], which further 
supports their topological nature. It should be noted that the quadratic f{R) model considered here also 
exhibits nonsingular cosmological solutions with a bounce [48-50] (see also [32,51] for related results). 
Whether such bouncing solutions can be related to the interior of the black hole solutions obtained here 
is an issue that will be explored elsewhere. 


5. Geodesics 


The tangent vector for null and time-like geodesics (/c = 0,1 respectively) parameterized by an affine 
parameter A can be written using Equation (9) as [28,52] 

2 1 \ 2 / , \ 2 


k = 


A{x) 


+ 


dx 


dX 


^ (^) -XT 


dip 


dX 


( 22 ) 


Ir \dXj ' A{x)fR 

where the angular variable 9 has been set to 7 r /2 without loss of generality. Given the invariance of the 
metric under time translations and p rotations, we have two conserved quantities, namely, E = and 
L = For time-like geodesics, the constants E and L represent the energy and angular momentum 
per unit mass, respectively. Using these conserved quantities, the relevant equation becomes 




(23) 


Instead of attempting to integrate this equation directly, it is convenient to replace the variable x by 

1 

r = TcZ using the relation = zf^, which leads to an equation that relates 2 ; and A = X/rp 


dX 

dz 


/I 


1 + 


2/fl 


(24) 


E2f2 


R 


A{^)fR 


k -f 


From the positivity of the argument in the square root of the denominator of this equation, as well 
as from the right-hand side of Equation (23), it is rather obvious that for fc = 1 and/or L 7 ^ 0, all such 
geodesics will bounce at a finite value 2; > 1, which corresponds to the location at which dx/dX = 0 
in Equation (23). Note, in this sense, that as 2; —>■ 1 the E^f"^ term goes to zero much faster than 
the other term in the square root. In fact, the leading order term of A{z)fR {k -f as 2: —)■ 1 is a 

positive constant, ’ which necessarily implies the vanishing of the square root at some 

z > 1 for arbitrary E. A similar behavior is also found in the standard Reissner-Nordstrom solution 
of GR, where time-like geodesics and null geodesics with nonzero angular momentum cannot reach the 
central singularity due to an infinite potential barrier (see [52] and [53] for details). Null geodesics, 
however, do reach the central singularity in the Reissner-Nordstrom case, and do it in a finite length 
of the affine parameter, implying in this way the existence of inextendible geodesics. As a result, the 
Reissner-Nordstrom space-time is regarded as singular. 

In our wormhole scenario, null radial geodesics satisfy the simpler equation 


±E 


dX 

dz 


1 + 


zfR,z 

2/h 


I 


(25) 
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which admits an exact analytical solution of the form 


±E\{z) 


yjl — z-2(“+l) 


+ 2Z2F1 


1 - 


2(a + l)’ 2(a + l) 




-2(a+l) 


(26) 


The asymptotie behaviors of this solution are the following: when z ^ 1, then ±EX ^ z, and when 
z ^ 1 we find ±EX —-^=1-^= = — jlj (see Figure 2). The divergenee of the affine parameter on 
both limits, when z 00 and when z 1, implies that null radial geodesies are eomplete. Note that in 
the limit z 1, light rays take an infinite (affine) time to reach the wormhole throat, where the eurvature 
divergenee is loeated. As a result, this spaee-time must be regarded as nonsingular, despite the existenee 
of eurvature divergenees at the wormhole throat. 


A(z) 



Figure 2. Representation of outgoing (green) and ingoing (red) null radial geodesies with 
E = 1 and a = 4/5. Similar results are found for any other value 0 < a < 1. 

6. Summary and Conclusions 

In this work, we have studied statie, spherieally symmetrie geometries generated by a partieular type 
of f(E) theory of gravity formulated a la Palatini eoupled to a one-parameter family of anisotropie 
fluids, whieh have also been put into eorrespondence with speeifie types of nonlinear theories of 
electrodynamies. We have shown that the resulting solutions represent wormholes and that their geodesie 
strueture is rather different from that diseussed in [18,52] for the Born-Infeld and Rieei-squared gravity 
models. In that ease, null rays and time-like observers (in eertain eonfigurations) were able to go through 
the wormhole. In the ease eonsidered here, however, neither light rays nor time-like observers ean ever 
reaeh the wormhole throat. In faet, we have shown that for null geodesies the wormhole throat lies at 
an infinite affine distanee. This point is very important beeause eurvature sealars generieally diverge 
at the throat. Given that this region lies beyond the reaeh of physieal observers and light rays, the 
spaee-time must be regarded as nonsingular despite the existenee of eurvature divergenees. This provides 
a different way to get rid of spaee-time singularities whieh eomplements the results presented in [18,52] 
and confirms the versatility of Palatini theories of gravity to get rid of singularities in simple and elean 
seenarios. This is in sharp eontrast with the diffieulties found in the standard metrie approaeh, where 
the higher-order eharaeter of the field equations largely prevents the eonstruetion of analytieal solutions, 
unless additional eonstraints are imposed upon the theory and/or the solutions. Following the analysis 
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carried out in [18] it would be interesting to explore the propagation of waves in this f\R) seenario. 
Sueh study will be earried out elsewhere. 

From our results, it follows that the emergenee of the wormhole is a strongly non-perturbative 
phenomenon. In faet, a glanee at Figure 2 puts forward that the linear behavior A (2;) ~ 2; of GR for 
null radial geodesies is valid in the f{R) case from infinity down to the seale 2; ~ 2, with the wormhole 
loeated at 2: = 1. This means that there are no loeal experiments that can be performed away from the 
wormhole that may be used to infer its existenee, i.e., there is no way to distinguish between this f{R) 
theory and GR by looking at perturbative eorreetions. One must get right at the throat of the wormhole 
to pereeive its existenee. An important lesson we extraet from this result is that GR ean be trusted at high 
and very high energies beeause there are examples of theories (like the one eonsidered here) in whieh its 
shorteomings ean be resolved by sudden non-perturbative effeets. 
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